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ABSTRACT 

We describe a project (transitsearch.org) currently attempting to discover transiting 
intermediate-period planets orbiting bright parent stars, and we simulate that project's perfor- 
mance. The discovery of such a transit would be an important astronomical advance, bridging the 
critical gap in understanding between HD 209458b and Jupiter. However, the task is made diffi- 
cult by intrinsically low transit probabilities and small transit duty cycles. This project's efficient 
and economical strategy is to photometrically monitor stars that are known (from radial veloc- 
ity surveys) to bear planets, using a network of widely-spaced observers with small telescopes. 
These observers, each individually capable of precision (1%) differential photometry, monitor 
candidates during the time windows in which the radial velocity solution predicts a transit if 
the orbital inclination is close to 90°. We use Monte Carlo techniques to simulate the perfor- 
mance of this network, performing simulations with different configurations of observers in order 
to optimize coordination of an actual campaign. Our results indicate that transitsearch.org can 
reliably rule out or detect planetary transits within the current catalog of known planet-bearing 
stars. A distributed network of skilled amateur astronomers and small college observatories is a 
cost-effective method for discovering the small number of transiting planets with periods in the 
range lOd < P < 200d that orbit bright (V < 11) stars. 

Subject headings: planetary systems — solar neighborhood 



1. Introduction 

Over the past seven years, Doppler radial veloc- 
ity (RV) measurements have led to the discovery 
of over one hundred planets within a sample of 
several thousand bright, nearby Sun- like stars. As 
the catalog of worlds continues to grow, our view 
of extrasolar planets is shifting from an anecdotal 
collection of individual systems, e.g. 51 Pegasi, v 
Andromedae, or 47 Ursae Majoris, to a more com- 
plete statistical census, in which categories and 
populations of planets can be clearly delineated 



(Marcy, Cochran, & Mayor 2000). 1 Yet the plan- 
etary systems from which we can learn the most 
— those that transit — remain anecdotal at best. 

For each system, there is a chance that the 
planet will periodically transit the surface of the 
star as seen from Earth. An eclipsing Jupiter-mass 
planet on a 3-day orbit produces a periodic ~1.5% 
dimming of the parent star that lasts for about 
3 hours. At present (August 2003) only a sin- 

An up-to-date version of the planetary census can be found 
at http:/ /cfa-www. harvard.edu/planets/ 



gle transiting planet — HD 209458b, P = 3.525d 
- has been studied in detail (Charbonneau et 
al. 2000; Henry et al. 2000), while a second ob- 
ject (OGLE TM-56-b, Konacki et al. (2003)) has 
been recently announced, but not studied as ex- 
tensively. HD 209458b has provided a scientific 
bonanza, including direct and accurate measure- 
ments of the planet's radius (1.35 ± 0.06i?j up ; 
Brown et al. (2001)), mass (0.69 ± 0.05M Jup ; 
Mazeh et al. (2000)), density, and even sodium 
in its atmosphere and hydrogen in its exosphere 
(Charbonneau et al. 2002; Vidal-Madjar et al. 
2003). 

The excitement generated by HD 209458b has 
led to a major push by the community to find ad- 
ditional transiting planets. A website maintained 
by Keith Home 2 lists, along with the project de- 
scribed in this paper, an additional twenty-four 
ground-based collaborations that are engaged in 
various efforts to discover planetary transits. In 
total, these surveys yield a reported capacity for 
discovering 148 planets per month. Despite this 
activity, however, an important corner of parame- 
ter space receives extremely little coverage; there 
is currently no other organized effort to detect 
intermediate-period planets which transit bright 
(V < 11 parent stars). We describe a strategy for 
detecting such transits, which we have adopted for 
the transitsearch.org collaboration, and we show 
Monte Carlo simulations that demonstrate the 
project's feasibility. 

Our basic approach is to harness a network 
of small independent telescopes to obtain multi- 
ple differential-photometric time-series of known 
planet-bearing stars during the well-defined time 
windows in which transits are predicted to occur. 
If several independent observers simultaneously 
measure a characteristic diminution or brighten- 
ing at the predicted times of ingress or egress, then 
there is strong evidence that the star is exhibiting 
a transit, and follow-up confirmation can then be 
obtained at the time of the next predicted tran- 
sit. The observational campaign is coordinated 
through a website: www.transitsearch.org. 

2 http://star- www.st-and.ac.uk/~kdhl/transits/tablc. html 



2. Scientific Motivation 

Any transits which our network uncovers will 
occur for planets which occult bright (V < 11) 
stars. This is an advantage. Such stars are pre- 
cisely those for which the RV method can provide 
accurate orbital parameters and accurate values 
for Msin(i), both of which are required to use- 
fully characterize the planetary properties. Fur- 
thermore, a bright parent star facilitates accurate 
photometry. The exquisite precision (1.1 x 10~ 4 ) 
HST light curves which have been produced for 
HD 209458b (Brown et al. 2001) depend on the 
V = 7.64 magnitude of the parent star. Brown 
et al. (2001) report that in order to obtain opti- 
mal photon noise-limited precision with HST for 
HD 209458b, photometric measurements of 80s 
duration (60s integration plus 20s CCD readout) 
were required. The critical ingress and egress pe- 
riods were thus time-resolved into approximately 
20 samples each. A V ~ 9 star, which produces 
^6 times fewer photons, would require 6-minute 
cadencing to obtain the same photometric preci- 
sion, and the periods of ingress and egress would 
be resolved into only a few time intervals. For con- 
siderably dimmer stars (V ~ 14, say) photometric 
precision will necessarily be compromised. 

The transitsearch.org collaboration is geared to 
survey planets with lOd < P < 200d. This sen- 
sitivity to longer-period transits occurs because 
we can narrow our observations to specific pre- 
dicted time windows. The detection strategy thus 
involves no data folding, and does not demand sta- 
ble photometry over multiple nights or seasons. 

Why would an intermediate-period transiting 
planet be of interest? Although the measured 
(1.35i?j up ) radius of HD 209458b is broadly con- 
sistent with its being a gas-giant composed pri- 
marily of hydrogen (Guillot et al. 1996; Burrows 
et al. 2000), recent work by Guillot & Showman 
(2002); Bodenheimer, Laughlin, & Lin (2003); 
Baraffe et al. (2003) all suggest that our under- 
standing of irradiated giant planets is incomplete. 
These three studies agree that standard evolu- 
tionary models can recover the observed radius of 
HD 209458b only if the deep atmosphere is unreal- 
istically hot. The recent studies incorporate real- 
istic atmospheric temperature profiles; the no-core 
models of HD 209458b have a radius of ^l.li?j up , 
which is much too small. 
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Three resolutions to this problem have been 
suggested. Bodenheimer, Laughlin, & Lin (2003) 
show that HD 209458b might be receiving interior 
tidal heating through ongoing orbital circulariza- 
tion resulting from perturbations due to a second 
planetary companion, whereas Guillot & Show- 
man (2002) propose that strong insolation-driven 
weather patterns on the planet are leading to con- 
version of kinetic wind energy into thermal energy 
at pressures of tens of bars. Burrows, Sudarsky 
& Hubbard (2003) argue that the size discrepancy 
stems largely from improper interpretation of the 
transit radius, and that the measured radius of 
HD 209458b in fact lies much higher up in the 
planetary atmosphere than is generally assumed. 

In any event, an accurate size and mass deter- 
mination for an intermediate-period planet will be 
of great help in resolving the observed size discrep- 
ancy for HD 209458b. A planet with intermediate 
period cannot not have significant internal tidal 
dissipation, but would still be receiving a mod- 
est amount of kinetic heating from the mechanism 
suggested by Guillot & Showman (2002) . Further- 
more, the Burrows, Sudarsky & Hubbard (2003) 
theory for HD 209458b is readily extended to pre- 
dict effective transit radii at different planetary 
masses and temperatures; the discovery of an in- 
termediate period transiting planet would provide 
a useful test of such predictions. 

Intermediate-period planets are also interest- 
ing because they can harbor dynamically sta- 
ble large satellites. Tidal interactions likely re- 
moved any satellites larger than R = 70km orbit- 
ing HD 209458b. Mars-mass moons, however, can 
last for 5 Gyr in the Hill Sphere of a lMj up planet 
orbiting a IMq star in a 27 day (0.18 Au) orbit, 
whereas in a 54 day (0.28 au) orbit, Earth-mass 
moons are dynamically stable (Barnes & O'Brien 
2002). Brown ct al. (2001) report that with HST, 
detections of satellites as small as 1R® are feasible. 
Therefore, an intermediate-period planet found by 
our survey could be followed up to search for large 
moons, and, additionally, planetary rings. Prior to 
space-based missions such as KEPLER (Borucki 
et al. 2003), the detection of a large moon orbit- 
ing an intermediate-period transiting planet is the 
best prospect for finding a habitable world. 



3. How Many Planets Transit Bright 
Stars? 

Regardless of scientific benefit, our survey can 
be successful only if there are additional transit- 
ing planets to be found orbiting bright stars, and 
only if the telescope network is sensitive and re- 
sponsive enough to definitively confirm or rule out 
the occurrence of transits for individual stars. 

The a priori probability that a planet transits 
its parent star as seen from the line of sight to 
Earth is given by 



= 0.0045 



lAU R+ — Rp 1 + ecos(-7r/2 - vo) 



Rp, 



1 -e 2 



(1) 

where a is the semi- major axis of the orbit, R± 
is the radius of the star, R p is the radius of the 
planet, e is the orbital eccentricity, and w is the 
argument of periastron referenced to the plane of 
the sky. Using the parameters of the current radial 
velocity planet catalog 3 , we find that among the 
17 Doppler wobble planets with periods P < lOd, 
there are (n) = 1.75 expected transits, and indeed, 
within this group, a transiting case (HD 209458b) 
is known. Sixteen of the planets with P < lOd 
have reported non-detections (although in some 
cases unpublished and unverified). These non- 
detections include HD 68988b, HD 168743b, and 
HD 217107b, which can be ruled out on the ba- 
sis of observations made with the transitsearch.org 
network. Only one P < lOd planet, HD 162020b 
(P = 8.428d), has, to our knowledge, not yet been 
checked for transits. 

Among the aggregate of 27 planets having pe- 
riods in the range lOd < P < 200d, the expected 
number of transiting planets is (n) = 0.72. Al- 
most none of the parent stars in this group, how- 
ever, have yet been monitored for transits, due 
to low individual transit probabilities and increas- 
ingly uncertain transit ephemerides. The main se- 
quence stars harboring intermediate-period plan- 
ets therefore represent the primary targets for our 
network. We also note that among the 67 known 
planets with P > 200d, one expects (n) w 0.6 ad- 
ditional transiting cases. However, as the plan- 
etary period becomes longer, follow-up becomes 
increasingly difficult due to uncertainties in the 
transit times and long intervals between occulta- 



3 see, e.g., http://www.transitsearch.org/stardatabasc/index.htm 
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tions. 

The transit probability for any given planet is 
not a strictly declining function of semi-major 
axis. For example, the highest a priori tran- 
sit probability for any known planet belongs not 
to one of the short-period hot Jupiters (which 
tend to average ^transit ~ 12%), but rather to the 
P = 550d planet orbiting i Draconis. In this sys- 
tem, a large planetary semi-major axis (1.34 Au) 
is more than offset by the 12.8i?0 stellar ra- 
dius (Allende Prieto & Lambert 1999), and fa- 
vorable orbital geometry (e=0.7, w ~ 94°; Frink 
ct al. (2002)) which lead to transit = 15.4%, with 
the next predicted transit occuring on April 4, 
2004. An extreme case such as this leads to 
a transit depth which is hard to detect from 
the ground (and impossible for our network) 
but there arc other intermediate-period planets 
which have surprisingly high transit probabilities 
(e.g. HD 38529b: P = 14.5d, transit = 13.7%, or 
HD 74156b: P = 51.6d, P transit = 4.3%). 

In addition to the current census, more plan- 
ets with periods suitable for transitsearch.org will 
emerge if RV surveys expand their samples. Cur- 
rently, within the lOd < P < 200d range, there 
are five known planets orbiting stars with V < 6, 
seven orbiting stars with 6 < V < 7, six orbit- 
ing stars with 7 < V < 8, and two planets each 
in the 8 < V < 9, and 9 < V < 10 ranges. If 
we assume that every available chromospheri- 
cally quiet main sequence dwarf with V < 6 has 
been adequately surveyed for P < 200d planets, 
and that each magnitude bin of unit width con- 
tains 1.8 x as many stars as available for bin 
(V — 1) (Cox 2000), then we expect that roughly 
9 + 16 + 29 + 52 + 94 = 200 detectable planets 
with P < 200d exist in orbit around stars with 
V < 11, indicating that close to 180 additional 
planets in this category can be detected using cur- 
rent RV techniques for bright stars. Statistically, 
this implies that 6 intermediate-period transiting 
planets orbit bright nearby stars. 4 The goal of the 
transitsearch.org network is to find one of these 
transits. 



4 A similar calculation shows that 6 additional short-period 
(P < lOd) transiting planets are likely to be orbiting bright 
stars. 



4. Transit Detection With Small Tele- 
scopes 

In the past several years, a number of amateur 
astronomers have detected the HD 209458b tran- 
sits and have shown that the ~ 1% diminution 
produced by a transiting Jovian planet is read- 
ily observable via differential photometry obtained 
with small (8-10 inch aperture) telescopes fitted 
with commercial-grade CCD detectors. A report 
of one of these observations (Arto Oksanen, Sky 
& Telescope, January 2001) raised a provocative 
question: Is it a realistic possibility for a network 
of small-college observatories and highly experi- 
enced amateurs to discover a new transiting sys- 
tem? If so, many small telescopes can be organized 
to maintain a time-intensive volunteer-based tran- 
sit survey of known planet-bearing stars during 
predicted transit epochs. 

In order to investigate the viability of detecting 
transits using low-cost equipment and software, 
we designed and documented an end-to-end pro- 
cedure which allowed us to observe an HD 209458b 
transit. Our demonstration observatory consists 
of a Meade LX-200 8-inch f/10 telescope fitted 
with a Santa Barbara Instruments Group ST7E 
765x510 pixel CCD. Pointing/imaging/guiding, 
standard image reductions, and aperture photom- 
etry are accomplished with a laptop computer 
running The Sky, CCDSoft, and MIRA AP 6.0 
software, respectively. These tools are all well- 
documented, reliable, relatively inexpensive, and 
familiar to amateur astronomers. 

With a focal reducer, the CCD image of a tar- 
get region covers 36' x 24', which is generally large 
enough to admit several V = 9 — 11 comparison 
stars for differential photometry In the case of 
HD 209458, a field star HIP 108793 (V = 8.33) 
is situated 12' away. We acquire the target field 
prior to the predicted start of ingress, and obtain 
successive 2s CCD exposures at a cadence of 35s 
per frame. The short 2s exposures are used to 
avoid pixel saturation by HD 209458. The small 
overall duty cycle is caused by the need to acquire 
auto-guiding images and to read out the CCD. Se- 
quences of 20 exposures are averaged together to 
produce composite measurements of the bright- 
ness of the stars in the field within 12 minute 
intervals. Using standard aperture photometry 
techniques, photoelectron counts from the target 
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Fig. 1. — Detection of the planet transiting 
HD 209458, using the portable observatory de- 
scribed in the text. The data were taken from 
Fremont, CA on the night of 19/20 October 2001. 
The solid curve is the model of Brown et al. (2001). 



star are compared to counts from the compari- 
son star(s) in the field. A transit manifests it- 
self by the characteristic changes in the brightness 
of the target star during the predicted times of 
ingress and/or egress. This phenomenon is shown 
in Figure 1, for HD 209458b during the transit of 
10/19/01. Photometric errors for each 12 minute 
bin are of order 0.003 magnitudes, and are domi- 
nated by atmospheric scintillation. We note that 
simple improvements, such as the use of a broad 
band, neutral density, or spot filter (Castellano 
2000) to increase open-shutter time or reduce the 
difference in brightness between the target and a 
typical comparison star, could considerably im- 
prove precision. 

For our purposes here, this observation tells us 
two things. First, transiting planets are readily 
detected with standard amateur-oriented equip- 
ment. Second, we can assume that many observers 
worldwide will have similar observational configu- 
rations, and will be capable of obtaining differen- 
tial photometry of comparable precision. 

5. Monte Carlo Simulation 

In order to evaluate the viability of a collaboration- 
based transit survey of known planet-bearing 
stars, we have performed a Monte Carlo study 
which models realistic incarnations of the transit- 
search. org network. This simulation is written in 



IDL 5 and makes heavy use of the IDL Astronomy 
User's Library 6 . 

The simulation is initialized with the follow- 
ing inputs: a list of observers and a target list of 
known planet-bearing stars. Each observer has an 
associated location (latitude and longitude) and 
weather (average fraction of clear/cloudy nights 
per year). Each target has an associated position 
(RA and Dec), period P, estimated transit prob- 
ability ^transit calculated from equation 1, and a 
transit ephemeris. The actual ephemeris for any 
given system can be obtained from fits to exist- 
ing RV data, but for simplicity in this simulation, 
transit ephemerides are generated randomly in- 
stead. In cases where a single star hosts multiple 
planets, it is listed in the target list with multiple 
entries. Once observers and targets have been set 
up, several record-keeping logs are also initialized. 

The first Monte Carlo step is to assign which 
targets will host real transits in the simulation. 
The program calculates a true/false condition for 
each target based on its ^transit- The simula- 
tion then enters its main loop which proceeds 
through an observing campaign night-by- night; 
within each night there is nested a loop which pro- 
ceeds through the observer list one-by-onc. That 
inner loop over observers proceeds as follows. 

Before assigning a target to an observer, the 
simulation first must determine if the weather 
is favorable for the night. While season-based 
weather patterns have not been figured into the 
model, the fraction of clear and cloudy nights at 
each observing location is known 7 . Using these 
probabilities, the night's weather for each location 
is determined in Monte Carlo fashion. If clear or 
cloudy, all observers common to the location will 
be affected identically. If the weather is deter- 
mined to be partly cloudy, each observer must be 
dealt with independently, allowing for the possi- 
bility that some observers in a given location will 
be able to observe while others are not. 

For each observer, the Julian date (JD) of sun- 
rise and sunset at the observer's location are cal- 
culated for the current date in the campaign. The 
JD of sunset, the position of each target, and the 

5 http://www.rsinc.com 
b http:/ /idlastro. gsfc.nasa.gov 

7 For locations in the US, such data are available from 
http://www.ncdc.noaa.gov/oa/climatc/online/ccd/cldy.html 
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observer's location are used to calculate airmasses 
for each target in the target list. Any targets that 
pass the airmass cutoff (2.5 in our simulations) at 
sunset are then checked to ensure they will pass 
the airmass limit for at least four hours. Thus a 
night of data will consist of at least four hours of 
time-series photometry, and we allow it to be as 
long as nine hours if the target is up (and the Sun 
is down). 

With the narrowed list of targets that are up, 
the simulation next checks to see if any targets are 
near transit. We assume that transit ephemerides 
are accurate to within ~5% of an orbit. If the ob- 
server's night overlaps at all with this margin of 
error for a target that is up, then that observer 
is assigned to that target. The likelihood of two 
observable stars being at transit in a given night 
is low, but if such a case does occur, the observer 
is assigned to the higher transit probability. Al- 
though there will be fewer opportunities to observe 
the longer-period planets, we feel it is justified to 
concentrate resources on targets where the proba- 
bility of successfully observing a transit is higher. 

Thus, only if an observer has a target near- 
transit and favorable weather, the simulation 
generates the night's photometry via Monte 
Carlo. Photometry is generated with an arbi- 
trary zeropoint and Gaussian noise. Addition- 
ally, for those targets that the simulation has 
randomly designated "real" transits, a simple 
linear ingress/egress and transit depth based on 
HD 209458b are input into the photometry. We 
scale the ingress/egress times with the period of 
the planet in order to simulate long-period tran- 
sits. The transit depth and photometry noise am- 
plitude are always set to match our template data 
(plotted in Figure 1) from the system described in 
§4- 

An observer's photometry for a night is very 
simply analyzed by calculating the Spearman's 
rank correlation for the data. Basically, this deter- 
mines whether a linear trend has been detected be- 
tween the beginning and end of a night's data and 
returns a confidence-level for such a trend (Press et 
al. 1992). This simple analysis fits the overwhelm- 
ing majority of cases where portions of ingress or 
egress have been observed, but will fail in the ex- 
ceedingly rare case where the transit is perfectly 
centered in the night's time window. 

The process of the preceding five paragraphs is 



repeated for each observer in the observer list, and 
this constitutes one night of the campaign (the 
inner loop). The simulation then increments its 
internal "calendar" by one day and proceeds again 
— this repeats until the campaign ends (the outer 
loop). 

For each target, the simulation keeps track of 
the number of no-correlations, 2-cr correlations, 
and 4-cr correlations observers have seen. Any 2-cr 
correlations are used to give a target a "free pass" , 
allowing it to stay on the target list until more 
definite observations can be made. However, since 
the goal of any campaign is to observe as much of 
the target list as possible, targets must be elim- 
inated from the list. Thus, a limit is set on the 
number of times a 4-cr correlation is shown before 
a star is dropped. This limit is generally low, for 
two reasons: in a real run, someone with access 
to a professional observatory will begin follow-up 
work, and additionally, in multiple realizations the 
4-cr correlation was associated with false positives 
in less than 1% of all cases. Similarly, a limit is 
set on the number of times a target can produce a 
non-detection (no-correlation) before it is dropped 
from the list. This limit is generally fairly high, to 
avoid dropping an actually-transiting target sim- 
ply because the transit occurred outside an observ- 
ing window. Both drop limits scale with target 
planet period, as the 5% accuracy of ephemerides 
leads to exceedingly long observing windows for 
targets with periods as short as 100 days, leading 
necessarily to a higher number of non-detections. 

In addition to this scorecard for each target, 
the simulation also generates many other record- 
keeping files, such as a log of the weather and 
observations of each observer for each night, and 
every photometry file that is generated. Certain 
variables of interest are also tracked, such as the 
summed transit expectation value from the re- 
maining target list, updated whenever a target has 
been successfully eliminated from further observ- 
ing. 

6. Results 

We have modeled several different configura- 
tions of observers, and a number of target lists. 
Addressed in this paper are five observer scenar- 
ios: a lone, dedicated observer at Mt. Hamil- 
ton near San Jose, CA; ten observers located at 
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Mt. Hamilton; ten observers distributed across 
the U.S.; twenty observers split between San Jose 
and Sydney, Australia; and twenty observers dis- 
tributed between eight worldwide locations. For 
comparison, a second run of a single observer has 
been made, but using a noise amplitude one-fifth 
that of the other runs, to represent an astronomer 
with access to a professional observatory. In each 
case, simulations are run for three target lists 
drawn from the pool of planet-bearing stars. The 
selection is based on planetary period, with max- 
imum period cutoffs of 1000 days (essentially the 
complete listing of extrasolar planets), 365 days, 
and 100 days. In each case, 7 days is always the 
minimum period cutoff. 

During trial runs, the 4-cr correlation and no- 
correlation limits were adjusted to minimize the 
number of false positives/negatives, while not 
making runs overly long. The values we used for 
these were such that there were no false positives 
among the amateurs, and two for the professional 
observer, for which the drop limits were lowered 
to match the reduced noise factor. While several 
actual transits were missed, these were due to in- 
completeness at the end of the run, not due to 
observers incorrectly ruling them out. 

For each pair of observer list and target list, the 
simulation has been run 20 times to provide ade- 
quate statistics, and was analyzed to find a num- 
ber of quantities, most importantly list complete- 
ness. Though it is possible to run the code until 
every target has been adequately observed, this 
usually leads to prohibitively long runs. Often, the 
end of a run is dominated by long period planets 
which provide few opportunities for observing, and 
depending on observer locations, may not ever be 
observed. Rather than running to absolute com- 
pletion, we chose a fixed length for the run. In 
each case presented here, runs start near the end 
of summer, 2003, and end in January 2008. This 
generally means that the observers will not have 
completed the entire target list, but still provides 
a good demonstration of the differences between 
observer configurations. 

The final list completeness is tallied in Table 1. 
Additionally, plots of targets remaining to be clas- 
sified vs. time are shown in Figure 2, where 
we have recorded the JD on which targets were 
dropped and the number of targets remaining af- 
terward. A smoothed average curve is superim- 



posed upon the target vs. time data from all 20 
runs. Number of observers increases from left to 
right, as does the average physical separation be- 
tween observers. The size of the target list in- 
creases downwards. The red curve is the test case 
for our professional astronomer. 

In all cases, the curves show a characteristic de- 
lay time, in which observers begin to classify tar- 
gets, but have not had enough nights to remove 
targets from consideration. At the ~2 year mark, 
targets begin to become saturated, and the list 
of remaining targets shortens rapidly. Eventually, 
most targets are removed, and the curve flattens 
out as only the most difficult targets remain. As 
one would expect, having more observers increases 
the number of targets it is possible to cover by 
the end of the run, and additionally, it is easily 
seen in the 100 day list and 365 day lists that 
observers more evenly spread in longitude do a 
slightly better job completing the target list. Also, 
note that while a lone amateur observer virtually 
never amasses enough data to drop a target, even 
an observer with full access to a professional-grade 
telescope can do no better than 10 observers clus- 
tered around the same location. 

The differences are better illustrated when we 
plot not the number of targets remaining, but 
the summed transit expectation value remaining 
on the target list (see Figure 3). It is easier to 
see the effect of longitudinal spread among ob- 
servers. The second and fourth columns repre- 
sent observers concentrated in one and two loca- 
tions, respectively, and show almost three years 
before significant target completeness begins to 
show. By contrast, the third and fifth columns 
represent observers located over a larger spread in 
longitude, and they begin to make progress fully 
a year before their counterparts. Globally, the 
optimal configuration seems to be a large num- 
ber of observers with maximal spread in longi- 
tude/latitude. From our data, we show that this 
configuration of observers can cover roughly 40% 
more of the probability-weighted target list than 
even the professional observer. 

7. Conclusion 

These Monte Carlo simulations are a feasibility 
study that demonstrates the efficacy of the tran- 
sitsearch.org project. They differ from reality in 
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Table 1 
Simulation Results 





Period upper limit for target list 


Observer Configuration 


lOOd 


365d 


lOOOd 


1 observer, San Jose 


0.5 ± 1.5 


0.4 ± 1.0 


0.5 ± 0.7 


1 professional observer, San Jose 


55 ± 7 


45 ± 4 


47 ± 4 


10 observers, San Jose 


62 ± 5 


52 ± 3 


41 ± 3 


10 observers, U.S.A. 


65 ± 6 


56 ± 3 


44 ± 3 


20 observers, San Jose & Sydney 


74 ± 6 


64 ± 4 


50 ± 4 


20 observers, worldwide 


95 ± 3 


76 ± 4 


56 ± 3 



Note. — Percent of target list completed by end of run for each ob- 
server configuration and target list. 



important ways. For instance, in the simulations 
a simple rank-correlation analysis is applied to in- 
dividual observers' data, for efficiency In reality, 
multiple observers' data will overlap, and arbitrar- 
ily sophisticated techniques (along with eyeballs) 
will be brought to bear on anything that appears 
interesting. A shortcoming of the simulations is 
that the photometric noise and transit signal are 
based on the data from our demonstration obser- 
vatory described in §4. Intermediate-period plan- 
ets will likely have smaller radii than HD 209458b, 
and hence the transit signal will not be as deep. 
However, many transitsearch.org "amateurs" ob- 
tain photometry comparable to data from our 
setup — and the simple improvements mentioned 
in §4 and Castellano (2000) will increase photo- 
metric precision. In addition, the longer timescale 
of intermediate-period transits will allow for more 
binning of the photometry, further increasing pre- 
cision. Lastly, the transit depth is also very sensi- 
tive to the stellar radius, which varies significantly 
in the target list, but has not been considered here. 
Nevertheless we feel that this feasibility study is at 
least a reasonable demonstration of our strategy. 

These simulations show that while a single ob- 
server campaign is capable of discovering transits, 
this observer will generally leave 30-50% of the sky 
uncovered. Not only can multiple observers better 
cover the sky, they can also cover it more quickly. 
Additionally, the data reveal the importance of 



having not only multiple observers in multiple lo- 
cations, but also in ensuring that the observers 
cover a wide range of longitudes in both hemi- 
spheres. Note, for instance, the difference in time 
to completion between the case where 10 observers 
are located in both San Jose and Sydney, and the 
case where 20 observers are scattered across nine 
worldwide locations. It is apparent that longitudi- 
nal coverage is important. One naturally expects 
that weather will be a key factor in determining 
time to completion, as it will most dramatically 
affect the length of a run that is confined to a sin- 
gle location. But the spread in longitude proves 
equally important, as one might guess from the 
process of viewing eclipses on Earth. Both timing 
and location are everything. 

Most current work on transits is divided into 
two categories, our Mount Hamilton case (the sin- 
gle dedicated observer), and studies like OGLE 
which rely on time sequenced, wide-field snapshots 
that detect possible transits. However, we have 
shown that a single observer is at a disadvantage, 
no matter how powerful the telescope, while wide- 
field surveys suffer from false positives associated 
with binary stars, and additionally provide poor 
targets for follow-up radial velocity work. In the 
end, even confining a search to the known extra- 
solar planets produces a long list of potential tar- 
gets that proves difficult to work through. We 
have shown that by handing the bulk of observ- 
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ing work to a dedicated team of observers with 
good longitudinal coverage, we may ensure that 
when a transit is expected to occur, there is al- 
ways someone watching, and that this team will 
prove competitive with any other transit search 
venture. 
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Fig. 2. — Simulated performance of sets of observers, showing number of targets remaining vs. time. The five 
columns of figures, from left to right, show survey results from (1) A single observer at Mt. Hamilton, CA, 
(2) ten observers on Mt. Hamilton, (3) ten observers spread across the continental US, (4) twenty observers 
divided between two locations, one in each hemisphere, and (5) twenty observers distributed world-wide. 
The three rows, from top to bottom, show results from period-limiting the target list at 100, 365, and 
1000 days, respectively. The red curve represents a single observer on Mt. Hamilton, but with access to a 
telescope providing photometry ~5 times as accurate. A point is generated for every date on which a target 
was dropped from the list; points from all 20 realizations are shown. The line is a smoothed average over 
realizations. 
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Fig. 3. — Simulated performance of sets of observers, showing remaining transit expectation value vs. time. 
The five columns and three rows of figures are as in Figure 2. 
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